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Abstract: The channeling of the ion recoiling after a collision with a WIMP changes the 
ionization signal in direct detection experiments, producing a larger signal scintillation or 
ionization than otherwise expected. We give estimates of the fraction of channeled recoiling 
ions in Csl crystals using analytic models produced since the 1960's and 70's to describe 
channeling and blocking effects. 
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1. Introduction 

Understanding the possibility of having channeled ion recoils in crystals is important in direct 
dark matter experiments measuring ionization or scintillation signals Q-Q. These experi- 
ments search for Weakly Interacting Massive Particles (WIMPs) composing the dark matter 
halo of our galaxy through the energy they deposit in collisions with nuclei within a crystal. 
Channeling occurs when the nucleus that recoils after being hit by a dark matter particle 
moves off in a direction close to a symmetry axis or symmetry plane of the crystal. Chan- 
neled ions suffer a series of small-angle scatterings with lattice nuclei that maintain them in 
the open "channels" between the rows or planes of lattice atoms. Thus they penetrate much 
further into the crystal than in other directions and give 100% of their energy to electrons 
(their quenching factor is Q = 1), producing more scintillation and ionization than they would 
produce otherwise. In scintillators like Nal (Tl) or Csl (Tl), channeling increases the observed 
scintillation light output corresponding to a particular recoil energy. In dark matter searches, 
Csl (Tl) crystals are used by the KIMS collaboration [|J] . In this paper we give upper bounds 
to the geometric channeling fraction of recoiling ions in Csl crystals. 

We proceed here in a similar manner as we did for Nal, a very similar crystal, in a previous 
paper || (we have also already considered channeling of recoiling ions in Si and Ge crystals in 
a different paper Q). We use a continuum classical analytic model of channeling developed 
in the 1960's and 70's, in particular by Lindhard pf1- ||lq| . In this model the discrete series of 
binary collisions of the propagating ion with atoms is replaced by a continuous interaction 
between the ion and uniformly charged strings or planes. The screened atomic Thomas-Fermi 
potential is averaged over a direction parallel to a row or a plane. This averaged potential, U 
or Up, is considered to be uniformly smeared along the row or plane of atoms, respectively, 
which is a good approximation if the propagating ion interacts with many lattice atoms in 
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the row or plane by a correlated series of many consecutive glancing collisions with lattice 
atoms. We are going to consider just one row or one plane, which simplifies the calculations 
and is correct except at the lowest energies we consider. 

In order for the scattering to happen at small enough angles so that channeling is main- 
tained, the propagating ion must not approach a string or plane closer than a critical distance 
r c or x c respectively (r is the transverse distance to the string and x is the distance perpen- 
dicular to the plane). The critical distance depends on the energy E of the ion and on the 
temperature of the crystal. Ions which start their motion close to the center of a channel (for 
example ions incident upon the crystal from outside), far from a row or plane, are channeled 
if the angle their trajectory makes with the row or plane is smaller than a critical angle, V'c 
that depends on the critical distance of approach r c or x c , and are not channeled otherwise. 
Nuclei ejected from their lattice site by WIMP collisions are initially part of a row or plane, 
so they start their motion from lattice sites or very close to them. This means that "blocking 
effects", namely large-angle interactions with the nuclei in the lattice sites directly in front 
of the recoiling nucleus site, are important. In fact, as argued originally by Lindhard jjf, in 
a perfect lattice and in the absence of energy-loss processes, the probability that a particle 
starting from a lattice site is channeled would be zero. This is what Lindhard called the "Rule 
of Reversibility." However, any departure of the actual lattice from a perfect lattice due to 
vibrations of the atoms, which are always present, violate the conditions of this argument 
and allow for some of the recoiling lattice nuclei to be channeled. 

There are several good analytic approximations of the screened Thomas-Fermi poten- 
tial and each leads to a different expression for the transverse continuum string and plane 
potentials, U(r) and U p (x) respectively. As in Ref ||] here we use Lindhard's expression, 
because it is the simplest and allows to find analytical expressions for the quantities we need. 
The transverse averaged continuum potential of a string as a function of r, relevant for axial 
channeling, was approximated by Lindhard as 

U{r) = E^l-\n[—^ + l\, (1.1) 

where C is a constant found experimentally to be C ~ v3 and ipf = 2Z\Z2e 2 / {Ed). 
Here Z\ and Z2 are the atomic numbers of the recoiling and lattice nuclei respectively, d 
is the spacing between atoms in the row, a is the Thomas-Fermi screening distance, a = 
0.4685A(Z 1 1/2 + Z2 /2 )~ 2/3 0, § and E = Mv 2 /2 is the kinetic energy of the propagating 



ion. In our case, E is the recoil energy imparted to the ion after a collision with a WIMP, 

E = ^-, (1.2) 

2M' v ; 

where q is the recoil momentum. The string of crystal atoms is at r = 0. The transverse 
averaged continuum potential of a plane of atoms, relevant for planar channeling, given by 
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Lindhard as a function of x is 



U p {x) = Ei? a 



x 



(1.3) 



where tp 2 = 27rnZiZ2e 2 a/E, n = Nd pc h is the average number of atoms per unit area, N is 
the atomic density and d pc h is the width of the planar channel. Also, the axial channel width 
d ac h is defined in terms of the interatomic distance d as ci ac h = 1/y/Nd, with N the atomic 
density. The plane is at x = 0. Examples of axial and planar continuum potentials for Cs 
ions propagating in the <100> axial and {100} planar channels of a Csl crystal are shown 
in Fig. ||. The potentials for Cs and I ions are practically identical, because Zq s ~ Zi (see 
Appendix A). 
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Figure 1: Continuum axial (black) and planar (green/gray) potentials for Cs ions, propagating in 
the <100> axial and {100} planar channels of a Csl crystal. The screening radius shown as a vertical 
line is ac s = 0.007785 nm (see App. A). 



For a "static lattice," that here means a perfect lattice in which all vibrations are ne- 



glected, the critical distances of approach r c and x c are given in the Eqs. 1A and L5 , expres- 
sions that were derived in Ref [9] . The critical distance for axial channeling is 



Tc{E) 




r— . 1 (l-3z/2) ^ 
v I +scos( 3 arccos (i - j 



(1.4) 



where z = 9E\d 2 / (8EC 2 a 2 ) . For planar channeling we will follow the procedure of defining a 
"fictitious string" introduced by Morgan and Van Vliet fl(], [l6| . They reduced the problem of 
scattering from a plane of atoms to the scattering of an equivalent row of atoms contained in 
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a strip of a certain width centered on the projection of the ion path onto the plane of atoms. 
Thus, 

x c {E) = f c (E), (1.5) 



where f c (E) is the critical distance obtained from Eq. |1.4| for the fictitious string. Along the 
fictitious row, the characteristic distance d between atoms needs to be estimated using data 
or simulations which are not available for a Csl crystal. As explained in Ref. ||], the choice 
of d equal to the average interdistance of atoms in the plane d p , i.e. d = d p , yields a lower 
bound on x c , which translates into an upper bound on the fraction of channeled recoils into 
planar channels. 

So far we have been considering static strings and planes, but the atoms in a crystal 
are actually vibrating with a characteristic (one dimensional rms) amplitude of vibration 
u\{T) which increases with the temperature T. In principle there are modifications to the 
continuum potentials due to thermal effects, but we take into account thermal effects in the 
crystal through a modification of the critical distances found originally by Morgan and Van 
Vliet and later by Hobler [|l(| to provide good agreement with simulations and data. For 
axial channels it consists of taking the temperature corrected critical distance r c (T) to be, 



r c (T) = Vr 2 c (E) + [ciMT)] 2 , (1-6) 
where the dimensionless factor c\ in different references is a number between 1 and 2 (see 



e.g. Eq. 2.32 of Ref. |LJ and Eq. 4.13 of the 1971 Ref. [jE]]). For planar channels, following 



Hobler [16| we use a similar equation 

x c (T) = v^PI+Mif, (i.7) 

where again C2 is a number between 1 and 2 (for example Barret [17] finds C2 = 1.6 at high 



energies, and Hobler [16| uses C2 = 2). We will mostly use c\ = C2 = 1 in the following, to try 
to produce upper bounds on the channeling fractions. 

In Appendix B it is shown that the variation of the lattice size with temperature, char- 
acterized by the variation of the lattice constant aj a t with temperature, has a negligible effect 
on the channeling fractions. This is why we ignore this effect (not only in this paper but also 
in our previous papers 0, §). 

We use the Debye model to account for the vibrations of the atoms in a crystal. The one 
dimensional rms vibration amplitude u\ of the atoms in a crystal in this model is ]8|, [15|] 



ui(T) = 12.1 A 



$(9/T) 1 



e/T 4 



+ 7 (MO) 



-i 



1/2 



(1.8) 



where M for a compound is the average atomic mass (in amu), i.e. for Csl, M = (Mq s +M\) /2, 
and T are the Debye temperature and the temperature of the crystal (in K), respectively, 
and $(x) = i /* tdb/{f* - 1) is the Debye function. With M Cs = 132.9 amu and Mi = 126.9 



amu, then M = 129.9 amu. We take the Debye temperature of Csl to be = 125 K || 18|. 
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The one dimensional rms vibration amplitude u\ in Csl is plotted in Fig. | function of 
the temperature T. The crystals in the KIMS experiment were kept at °C in 2007 @. 
Currently the operating temperature of the crystals is 20 °C |]20[| . The vibration amplitude 
is ui = 0.0141 nm at °C, and m = 0.0146 nm at 20 °C. 
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Figure 2: Plot of m(T) for Csl (Eq. 1.8 with M = (M Cs + Mj)/2) 



Using the temperature corrected critical distances of approach r c {T) and x c (T) (Eqs. [I 



and |l.7| ) or the static lattice critical distances r c and x c (Eqs. |L4j and |l.5|) , we obtain the cor- 
responding critical axial and planar channeling angles ip c (see Ref. || for details). Examples 
are shown in Figs. |3| to [|, for c\ = C2 = c and c = 1 or c = 2 as indicated. 

Fig. [| clearly shows how the critical distances and angles change with temperature for 
a Cs or I ion propagating in the <100> axial and {100} planar channels of a Csl crystal, 
with temperature effects computed with c\ = C2 = c = 1. At small energies the static critical 
distance of approach is much larger than the vibration amplitude, so temperature corrections 
are not important. As the energy increases, the static critical distance of approach decreases, 
and when it becomes negligible with respect to the vibration amplitude U\, the temperature 
corrected critical distance r c becomes equal to c\u\. In this case, since u\(T) increases with 
T, the critical distance r c ~ c\U\ becomes larger with T, and therefore the critical channeling 
angle becomes smaller. Notice that for the 100 channels, the widths of axial and planar 
channels are the same, d ac h = d pc h and r c = x c . 

Fig. [I] shows the same effects for the <111> axial channel. In this channel, the critical 
distance (the minimum distance to a string to maintain channeling) becomes larger than the 
radius of the channel at energies below a few keV, shown in the figures. This means that 
nowhere in the channel an ion can be far enough from the string of lattice atoms for channeling 
to take place (thus the critical channeling angle is zero) . The exact calculation of the energy 
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Figure 3: Static (green) and temperature corrected with ci = C2 = c = 1 (black) (a) critical distances 
of approach (and u\(T) in red) and (b) the corresponding critical channeling angles, as a function of 
the energy of propagating Cs or I ions (they are practically the same for both) in the < 100 > axial 
(black) and {100} planar (green) channels. Here d ac h = d pc h- 



< 1 1 1 > axial channel, c= 1 



< 1 1 1 > axial channel, c= 1 
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Figure 4: Same as Fig. but for the <111> axial channel. 



at which this happens would require considering the effect of more than a single row of atoms 
(which we do not do here) thus our results at these low energies are only approximate. Notice 
that for the 111 channels, the <111> axial and {111} planar channels do not have the same 
widths, d acn ^ d pc h> and we only show the critical distance and angles for the axial channel. 

Figs. ||(a) and ||(b) show the static and T-corrected critical distances and angles repec- 
tively at several temperatures for traveling Cs or I ions in the 100 axial and planar channels 
with c\ = C2 = c = 2. 

2. Channeling fractions 

In our model, a recoiling ion is channeled if the collision with a WIMP happened at a distance 
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<100> and {100} channels, c=2 
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Figure 5: Same as Fig. || but with c\ = = c = 2. 



large enough from the string or plane from which the ion was expelled. Namely, channeling 
happens if the initial position of the recoiling motion is r« > rj )TO i n or X{ > x^min for an axial 
or planar channel respectively. Here rj jm i n and a^min depend on the critical distance, r c or x c 
(and through it on the energy E and the temperature T) and on the angle 4> the initial ion's 
momentum makes with the string or plane. In Ref. H we obtained the following expressions 
for both distances: 

Ca 



1 exp(-2 sin 2 0/Vi) 



titan ( 



(2.1) 



1 



with r c given in Eq. 1.6 and 



ft 

2 



C 2 



1 2 



+ c 2 - ^ 



sin" 



+ c 2 



Xc_ 

a 



8111' 



mi 



dp tan ( 



(2.2) 



where x c is given in Eq. [1/7. 

We take the initial distance distribution of the colliding atom to be a Gaussian with a 
one dimensional dispersion ui, and to obtain the probability of channeling for each individual 
channel we integrate the Gaussian between the minimum initial distance and infinity (a good 
approximation to the radius of the channel; see Ref || for details). The dependence of 
these probabilities on the critical distances enter in the argument of an exponential or an erfc 
function. Thus any uncertainty in our modeling of the critical distances becomes exponentially 
enhanced in the channeling fraction. This is the major difficulty of the analytical approach 
we are following. 

In order to obtain the total geometric channeling fraction we need to sum over all the in- 
dividual channels we consider. Taking only the channels with lowest crystallographic indices, 
100, 110 and 111, we have a total of 26 axial and planar channels, as explained in Appendix 
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A of Ref. || (Csl and Nal have the same crystal structure). Here "geometric channeling 
fraction" refers to assuming that the distribution of recoil directions is isotropic. In reality, 
in a dark matter direct detection experiment, the distribution of recoil directions is expected 
to be peaked in the direction of the average WIMP flow. The integral over initial directions 
is computed using HEALPix [^IJ (see Appendix B of Ref. ||). 

Fig. |6] shows upper bounds to the channeling probability computed for each initial recoil 
direction direction q and plotted on a sphere using the HEALPix pixelization for (a) a E = 200 
keV and (b) a 1 MeV Cs ion at 20 °C with c\ = C2 = 1 assumed for the temperature effects. 
The red, pink, dark blue and light blue colors indicate a channeling probability of 1, 0.625, 
0.25 and zero, respectively. The results are practically identical for an I ion. 

Fig. [7| shows upper bounds to the channeling fractions of Cs recoils for individual channels, 
for T = 20° C and assuming c\ = C2 = 1. The black and green (or gray) lines correspond to 
single axial and planar channels respectively. The upper bounds of the channeling fractions 
of planar channels are more generous than those of axial channels because of our choice of x c 
in Eq. 1.5. This does not mean that planar channels are dominant in the actual channeling 
fractions. 





Figure 6: Upper bounds on the channeling probability of a Cs ion (for an I ion the figure would be 
practically identical ) as function of the initial recoil direction for a (a) 200 keV and (b) 1 MeV recoil 
energy at 20 °C (with c\ = c% = 1). The probability is computed for each direction and plotted on 
a sphere using the HEALPix pixelization. The red, pink, dark blue and light blue colors indicate a 
channeling probability of 1, 0.625, 0.25 and zero, respectively. 

Upper bounds to the geometric channeling fractions of Cs and I ions as function of 
the recoil energy are shown in Figs. || and |9| with thermal effects taken into account with 
ci = C2 = 1 and c\ = C2 = 2, respectively. 

Notice that we have not included here any dechanneling due to the presence of impurities 
in the crystal (such as Tl atoms), which would decrease the channeling fractions presented. 
We intend to return to this issue in a later paper. As we see in Fig. || and ||, when no 
dechanneling is taken into account, the channeling fraction increases with energy, reaches a 
maximum at a certain energy and then decreases as the energy increases. This maximum 
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Figure 7: Upper bounds on the channeling fractions of Cs recoils as a function of the recoil energy 
E when only one channel is open, for T = 293 K with temperature corrections included in the critical 
distances with the coefficients c\ = c% = 1. Black and green/gray lines correspond to axial and planar 
channels respectively. Solid, dashed, and dotted lines are for 100, 110, and 111 channels respectively. 
The corresponding figure for an I ion would be practically identical. 
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Figure 8: Upper bounds on the channeling fraction of Cs (solid lines) and I (dashed lines) recoils as 
a function of the recoil energy E for T = 600 °C (orange/medium gray), 293 K (green/light gray), 273 
K (black), and 77.2 K (blue/dark gray) in the approximation of c\ = c% = 1 without dechanneling. 



occurs because the critical distances decrease with the ion energy E, which makes channeling 
more probable, while the critical angles also decrease with E, which makes channeling less 
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Static lattice 




probable. At low E the critical distance effect dominates, and at large E the critical angle 
effect dominates. 

As shown in Fig. ||, the channeling fraction for Csl (Tl) is never larger than 5% at 293 K 
(with cl = c2 = 1) and the maximum fraction happens at around 1 MeV. This is comparable 
to the channeling fraction of Na (Tl) which is also never larger than 5%, but in the case of Na 
(Tl) the maximum happens at 100's of keV (see Fig. 14(a) of Ref. H). For Si and Ge (with 



cl = c2 = 1) the channeling fractions reach about 1% and the maximum happens at 100's 
of keV for Si and at around 1 MeV for Ge (see Fig. 18 of Ref. ||). Fig. |{| shows that the 
maximum channeling fraction for Cs or I recoils at 293 K would be below 0.5% if c\ = C2 = 2 
instead. However, since we do not know which are the correct values of the crucial parameters 
ci and C2 for Csl, we could ask ourselves how the upper bounds on channeling fractions would 
change if the values of these parameters would be smaller than 1 (please recall that for other 
materials and propagating ions the values of these parameters were found to be between 1 
and 2). The values of c\ and C2 cannot be smaller than zero, thus Fig. [H| shows our most 
generous upper bounds on the geometric channeling fraction, obtained by setting c\ = C2 = 0, 
namely by neglecting thermal vibrations of the lattice (which make the channeling fractions 
smaller as T increases) but including the thermal vibrations of the nucleus that is going to 
recoil (which make the channeling fraction larger as T increases). Although it is physically 
inconsistent to take only the temperature effects on the initial position but not on the lattice, 
we do it here because using ci = C2 = 0, namely a static lattice, provides an upper bound on 
the channeling probability with respect to that obtained using any other non-zero value of c\ 
or C2- Even in this case, the channeling fractions at 293 K cannot be larger than 10%. 

To conclude, let us remark that the analytical approach used here can successfully de- 
scribe qualitative features of the channeling and blocking effects, but should be complemented 
by data fitting of parameters and by simulations to obtain a good quantitative description 
too. Thus our results should in the last instance be checked by using some of the many 
sophisticated Monte Carlo simulation programs implementing the binary collision approach 
or mixed approaches. 
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A. Crystal structure and other data for Csl 

Csl is a diatomic compound that has two interpenetrating face-centered cubic (f.c.c.) lattice 
structures displaced by half of a lattice constant with 8 atoms per unit cell. The lattice 
constant of Csl crystal is a\ at = 0.45667 nm at room temperature (Table 3.4 in Appleton and 



Foti |L5|]). The temperature dependence of ai a t is explained in Appendix B. 

The atomic mass and atomic number of Cs and I are Mq s = 132.9 amu, Mj = 126.9 amu, 
Zq s = 55 and Zj = 53. 

With respect to the Thomas-Fermi screening distance, for Cs recoils from a mixed row 
or plane we use the average 

acs = (acsCs + a Cs i)/2 = 0.007785 nm, (A.l) 
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where a Cs c s = 0.4685(Zc{ 2 + ^{ 2 )" 2/3 = 0.007761 nm and a Cs i = 0M85(Z^ + Z r 1/2 )- 2/3 = 
0.007809 nm correspond to a Cs scattering from a Cs and an I lattice atom, respectively. On 
the other hand, for Cs recoils from a pure row or plane we use acsCs because the row or plane 
from which the recoiling Cs ion was emitted contains only Cs atoms. Similarly, for I recoils 
from a mixed row or plane, we use 

a : = (a n + a Cs i)/2 = 0.007833 nm, (A. 2) 

where a n = 0.4685(Zj 1/2 + z\ /2 )- 2 l z = 0.007857 nm and acsl correspond to an I ion scattering 
from an I and a Cs lattice atom, respectively. For I recoils from a pure row or plane we use 
an , since the row or plane the recoiling ion is emitted from is made of I ions only. 

To compute the interatomic spacing d in axial directions and the interplanar spacing d p 
in planar directions, we have to multiply the lattice constant by the following coefficients ||: 

• Axis: < 100 >: 1/2 , < 110 >: l/y/2 , < 111 >: y/3/2 

• Plane: {100} : 1/2 , {110} : 1/2^2 , {111} : 1/2^3 

The Debye temperature of Csl is = 125 K, and the crystals in the KIMS experiment 



are currently at a temperature of 293 K [20|. 



B. Temperature dependence of lattice constant 

In general the lattice constant, ai at is temperature dependent. The change in a\ &t with 
temperature depends on the thermal expansion coefficient, (3 of a crystal. For Csl (Tl), 
(3 = 54 x 10~ 6 °C~ 1 . To find the change in the lattice constant at a temperature T and the 
lattice constant at 20 °C, we have 

[ai at (T) - a lat (20 °C)]/a lat (20 °C) = /3(T - 20 °C), (B.l) 

where ai a t(20 °C) = 0.45667 nm is the lattice constant at 20 °C for Csl (Tl). When T changes 



from 20 °C to 600 °C, the change in the lattice constant (using Eq. Rl) is only 3.1%. This 
change in a\ a t between 20 °C and 600 °C results in a negligible change in the channeling 
fractions. As an example the Cs channeling fractions with c\ = ci = 1 and c\ = C2 = 2 for 
the two choices of ai a t(20 °C) and ai at (600 °C) are shown in Fig. |i~T|(a), As the two curves 
are very similar, we use ai at (20 °C) for all three crystal temperatures in this paper. 



We can use Eq. B.l to find the temperature dependent lattice constant for Nal (Tl), Si 
and Ge crystals. For these three crystals, the coefficient of thermal expansion is /?Nai(Ti) = 
47.4 x 10~ 6 °C-\ /3 S i = 2.6 x 10~ 6 °C" 1 , and (3 Ge = 5.9 x 10~ 6 °C- 1 . When T changes from 
20 °C to 600 °C, the change in the lattice constant of Nal (Tl) is 2.75%. In Si and Ge, we can 
go to higher temperatures, and the maximum temperature that we considered in our previous 
paper on Si and Ge || was 900 °C. When T changes from 20 °C to 900 °C, the change in the 
lattice constant of Si and Ge is 0.23% and 0.52% respectively. 
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Csl Crystal, T=600 °C 



Nal Crystal, T=600 °C 
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Figure 11: Channeling fraction of (a) Cs ions propagating in a Csl crystal and (b) Na ions propagating 
in an Nal crystal as a function of the recoil energy E for T — 600 °C with ai a t(600 °C) (black) and 
ai a t(20 °C) (green/gray) for the two choices of ei = c 2 = 1 or 2. 



Si Crystal, T=900 °C Ge Crystal, T=900 °C 
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Figure 12: Channeling fraction of (a) Si ions propagating in a Si crystal and (b) Ge ions propagating 
in a Ge crystal as a function of the recoil energy E for T = 900 °C with ai a t(900 °C) (black) and 
aiat(20 °C) (green/gray) for the two choices of ci = c% = 1 or 2. 

The Na channeling fractions with c\ = c% = 1 and c\ = ci = 2 for the two choices of 
ai at (20 °C) and ai at (600 °C) are shown in Fig. |ll|(b) for an Nal crystal. Fig. |l2| shows the 
channeling fractions for Si and Ge for the two choices of ai a t(20 °C) and ai a t(900 °C). 

Clearly, the change in the curves is negligible. Thus we always used the value of ai at 
measured at 20 °C in this paper as well as in our previous papers || [|. 
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